function Wdraw=invWishartDraw(v,SInv)
% 
% [logPDFNoConstant,logPDF]=invWishartDraw(v,SInv) 
% 
% Generates a draw from the Inverse Wishart IW(v,S) 
% 
% *Note:* Input is Sinv=inv(S) 
%
%% Input 
% Sinv: [q,q] inv(S) Prior precision 
% v: [scalar] degrees of freedom
% 
% Alejandro Justiniano August 30 2013 (C) 
%==========================================================================
q=size(W,1); 
drMat=mvrndn(zeros(q,1),SInv,v); 
Wdraw=(drMat'*drMat)\eye(q); 
